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In  this  report,  the  propagation  errors  of  the  initial  con¬ 
dition  generator  are  studied  in  two  ways:  CO  by  comparing 
observational  data  on  one  pass  with  predictions  from  initial 
conditions  generated  on  a  previous  pass  of  the  same  satel¬ 
lite,  and  (2)  by  statistically  comparing  the  elements  gen¬ 
erated  by  different  passes  of  the  satellite.  These  two 
methods  are  found  to  agree.  Excluding  one  set  of  poor  data, 
it  is  found  that  the  propagation  error  has  a  position  stand¬ 
ard  deviation  of  about  ten  nautical  miles  perpendicular  to 
the  orbit  plane,  and  of  about  fifty  nautical  miles  in  this 

tion  of  about  seventy-five  to  one  hundred  seconds.  A  means 
of  reacquiring  the  satellite  on  the  basis  of  this  analysis 
is  presented  and  linearized  equations  for  the  antenna  pro¬ 
gram  are  listed. 
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FOREWORD 

This  Technical  Documentary  Report  was  prepared  under  Definitive 
Contract  AF04{695) -113  and  is  submitted  in  accordance  with  Exhibit  "A" 
of  that  contract  and  Paragraph  4.2.1  of  AFBM  Exhibit  58-1,  "Contractor 
Reports  Exhibit,"  dated  1  October  1959,  as  revised  and  amended. 

The  report  was  prepared  by  the  Philco  WDL  l^athemat ical  Analysis 
Department  in  fulfilling  the  requirements  of  Section  TI,  Tab  A,  Para¬ 
graph  1.2. 1.2  and  Section  II,  Tab  J,  Paragraph  13.1  of  AFSSD  Exhibit 
61-27A,  "Satellite  Control  Subsystem  Work  Statement,"  dated  15  February 
1962,  as  revised  and  amended. 

This  interim  report  was  prepared  to  partially  satisfy  the  require¬ 
ments  of  the  High  and  Low  Altitude  Tracking  Simulation  Studies. 
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AN  ACCWiACY  STUDY  OF  TiLt;  INITIAL  CONDITIONS  GENEiUTOU 
USING  VaiLOUT  DATA 


I.  INTliODUCTION 

The  "Initial  Conditions  Generator"  is  a  Fhilco  2000 
program  which  estimates  the  six  elliptic  elements  describing 
the  trajectory  of  an  artificial  satellite  from  a  set  of  observa¬ 
tions  of  the  satellite  on  a  single  pass  over  a  single  station. 

The  observations  may  be  either  Verlort  or  angle  plas  doji])ler. 

The  program  is  described  in  the  reference. 

The  estimated  elements  may  be  used  as  a  starting  point 
for  differential  correction.  In  this  case,  the  elements  will  be 
of  succificent  accuracy  if  the  differential  correction  converges. 
Since  the  differential  correction  usu.illy  converges  after  a  few 
iterations,  the  initial  conditions  generator  is,  indeed, 
adequate. 

However,  the  estimated  elements  may  also  be  used  to 
predict  an  ephemsris  dii  .itly.  This  elimination  of  the 
differential  correction  procedure  involves  a  significant  saving 
of  computer  time  at  the  exjiense  of  }>redlctioo  accuracy.  This 
note  describes  a  numerical  study  of  the  accuracy  of  predictions 
made  from  elements  estimated  on  the  basis  of  Verlort  tracking 
data. 

II.  DESCRIPTION  OF  METHOD: 

Several  passes  of  a  program  46l  and  a  Program  l62  satellite 
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over  various  Veriort  stations  have  been  used,  iiorae  characteris¬ 
tics  of  these  passes  arc  listed  in  Table  1.  Eac!:  of  these 
passes  was  run  tlirough  the  initial  conditions  generator  program 
to  obtain  elements.  Each  element  set  corresponding  to  the  Program 
461  satellite  was  then  run  through  the  *Generul  Predicitiona" 
program;  theProgram  l62  data  was  treated  in  a  different  manner 
described  in  section  V.  The  generul  x'fcdictions  program 
numerically  integrated  the  ecjuations  of  motion;  tested  whether 
the  satellite  was  in  the  field  of  view  of  one  or  more  observing 
stat  ons;  and,  if  so,  printed  the  range,  azimuth,  and  elevation 
of  the  satellite  at  each  time  step,  for  as  long  as  the 
satellite  remained  observable.  By  comparing  these  printed 
f;uantitie8  with  the  observations  actually  recorded  on  suc¬ 
cessive  ])US8es  of  the  satellite,  a  measure  of  the  propagated 
error  of  tlie  initial  conditions  gc.:erator  may  be  obtained. 

This  is  possible  because  the  propagation  error  is  much  larger 
than  either  the  integration  error  of  general  predictions,  or 
the  laudom  error  in  the  observations. 

Ill.  iVILD  POINT  itaiOVAL: 

Besides  the  occasional  wild  observations  usually  recorded 
by  electronic  sensors,  some  of  the  passes  used  iu  tliis  study 
coutain  caiihration  data.  On  t!i08e  passes  iu  fable  I,  for 
^Yhich  the  column  marked  "CiVL?"  is  checked,  the  measurements 
taken  in  the  course  of  calibrating  the  Veriort  are  trane- 
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milted  in  a  format  indiatinguishable  from  that  of  the  real 
obeorvations.  Since  these  calibration  meusurements  are  wildly 
inconsistent  with  the  orbit,  a  screening  procedure  is  necessary. 

Fortunately,  the  initial  condition  generator  includes  an 
optional  screening  procedure.  This  scheme  consists  of  cal¬ 
culating  the  average  linear  distance  between  the  observed 
points  and  the  plane  of  best  fit;  rejecting  those  points  which 
are  farther  from  the  plane  than  a  constant,  alpha,  times  this 
average  distance;  and  then  refitting  a  plane  to  the  remaining 
points.  After  some  experimentation,  it  was  found  that  this 
procedure  will  remove  both  the  calibration  and  the  isolated 
wild  points.  For  the  46l  data,  alpha  was  set  at  3>  and  for 
the  162  data  at  1.6.  In  both  cases,  the  rejection 
was  i  e]>eated  a  second  time.  The  reason  for  this  second  re¬ 
jection  is  that  the  average  distazice  calculated  the  first  time 
is  (juito  large,  so  only  the  most  horrible  points  are  rejected. 

IV.  i  iltlui)  KhUOK; 

It  is  a  well-known  but  rarely-published  fact  that  the 
bulk  of  tlje  propagation  ei'ror  arising  from  a  Bint,le  pass  orbit 
estimation  is  due  to  the  error  in  the  period.  That  is,  if 
the  period  is  estimated  exactly,  and  the  other  elements  are 
subject  to  small  errors;  the  propagated  error  will,  to  a 
first-order  theory,  be  a  periodic  function  of  time.  On  the 
other  hand,  a  small  error  in  the  period  will  produce  a  propaga- 
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tion  error  wlioye  magnitude  is  lou-hly  ^'roportional  to  prediution- 
time.  ,Ui  eijuivalont  interpretation  is  that  tlie  predicted 
position  and  velocity  co-ordinates  are  roughly  correct,  but 
that  the  associated  time  is  in  error  by  the  amount  ,  At. 


AP 


[1]  At-  f  (tp-  tj) 


where  P  is  the  period;  A  P,  the  period  error;  tp,  the  present 
time;  t  the  time  of  the  initial  conditions. 

^  i 

The  interpretation  just  given  to  the  inertial  co-ordinates 
cannot  be  directly  extended  to  the  radar  co-ordinates  -  - 
range,  azimuth,  and  elevation.  The  reason  is  that,  although  the 
predicted  position  at  time,  tp,  is  roughly  equal  to  the  true 
position  at  time,  tp  *  At;  the  position  of  the  observing  station 
is  different  at  these  two  times,  due  to  the  rotation  of  the 
earth.  Thus,  the  predicted  radar  co-ordinates  at  time,  tp , 
are  to  be  compared  with  the  following: 


[2J 


Sp  (Ap)  Sjj  (tp  +  At)  +  AS 
Ap  (tp)  (tp  +  At)  +  AA 

Ep  (tp)  (tp  +  At)  +  AE 
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where  a  'P'  subscript  iniiicuteu  pi  edicti  ii,  ana  aii  '  O' tj.iuscr ij-.t 
indie;. los  observation.  The  deltas  uay  be  a^.jiroximated  as 
liiiear  lunctions  of  it.  After  carrying  out  the  necessary 
differentiation 

[3]  is  At  R  c«8  E  sin  A  cos  L 

£1 

iA  At  (sin  L  -  cos  A  cos  L  sec  E  (sin  £  R/s)  ) 

AE  iJg  At  sin  A  cos  L  (l  +  R  sin  E/s), 

where  £1_  ia  the  rotation  rate  of  the  earth  (expi  essed  in 
£ 

radians  per  second  for  AS)  and  in  degrees  per  second  for  the 
angles);  L;  "the"  latitude  of  the  station;  and  it,  "the"  radius 
of  the  earth.  Then  equations  [2]  and  [3]  may  be  solved 
simultaneously  for  At.  Indeed  three  values  of  At  may  be  de¬ 
rived  and  averaged  at  each  prediction  time.  Finally,  from 
tj,  and  F,  AF  may  be  calculated. 

Equation  [3]  may  also  be  used  in  the  field  as  an  antenna 
dii  uctin^  pro^^nun  as  sut^^osted  in  the  conclusion. 

V.  DL.^CT  ZaN/iLYblb  OF  THE  ELEMENTS 

Instead  of  running  the  elements  estimated  by  the  initial 
conditions  generator,  through  g eneral  predictious,  it  is 
possible  to  analyze  them  directly.  Of  the  six  elements  -  one, 

'H  ,  is  a  rapidly  changing  function  of  time.  The  other  five 
elements,  suffer  only  perturbative  changes.  Thus,  if  the 
perturbations  are  sufficiently  small,  the  element  sets 
estimated  for  different  passes  of  the  same  satellite  may  be 
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treated  as  random  aaiiii)les.  In  particular,  an  estimate  oi‘  the 
i-eriod,  may  be  derived  from  the  other  elements  as  follows: 

0  2 

[4J  e  =  (r  +  W^) 

A  =|l{1  -  e2)j 
P  =  2  n 

wliere  ji  is  the  product  of  earth's  mass  and  the  universal  gravita¬ 
tion  constant.  I'or  the  46i  data,  the  estimates  of  'P  '  cal¬ 
culated  by  eijuation  E4]  ,  and  the  estimates  of  'ip'  calculated 
by  the  procedure  of  section  IV  have  been  obtained.  The 
quantity,  'P  -  iP'  was  then  computed  i or  each  pass.  This 
quantity  is  an  estimate  of  the  true  period;  its  value  should 
be  approximately  the  same  for  ouch  jiass.  In  table  II,  the 
values  of  'P 'iP'and  'P-AP'  are  listed,  along  with  their 
means  and  standard  deviations.  It  should  be  noted  that  the 
standard  deviation  of  P-AP  is  much  smaller  than  that  of 
the  other  two  -  -  thus,  the  procedure  is  self-cousistent. 

Because  ol  this  self-consistency,  it  was  deemed  unnecessary 
to  apply  the  general  prediction  pi’Ogram  to  the  162  data. 
The  btaiidard  deviation  of  the  jicriod  error  was  obtained  from 
the  disi^er.'ion  of  the  estimated  periods  calculated  by 
equation  [ 4  j  . 

In  tables  III  and  IV,  the  estimates  of  five  elements  -  - 
8emi-.Tjajor  axis,  A;  eccentricity,  right  ascension  of  ascending 
node,  ft  ;  inclination,  I;  and  j.eriod,  P  —  together  ith  the 
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meuiib  ;-au  .standard  deviations  are  listed  for  the  461  and 
162  satellites,  respectively.  'A',  and  'P'  are 

calculated  from  equations  [4]  ;  'I*  is  taken  directly  from  Lhe 
initial  conilitions  geuerutor is  also  taken  from  the 
initial  conditions  generator,  but  the  values  are  corrected  by 
a  linc.ir  function  of  time  to  remove  the  well-known  secular 
perturbation  caused  by  the  first  harmonic  of  the  eartli'a 
potential . 

VI;  CO.N'Cl  USIONri : 

A  glance  at  tables  II  and  III  shows  that  serial  number 
544  contains  much  jioorer  data  than  the  other  46l  passes. 

An  examination  of  the  data  listing  reveals  that  active 
tracking  was  lost  while  the  satellite  was  still  a, .preaching 
the  station.  Aside  from  this  peculiarity,  the  observations 
look  reasonably  smooth.  If  we  eliminate  this  one  bad  pass 
fronithc  analysis,  the  results  are  greatly  improved.  Table  V 
shows  the  means  and  standard  deviations  of  the  estimated 
46l  elements  as  calculated  from  the  four  good  passes  -  - 
551,  55-,  553 1  561.  It  sliould  be  noted  that  these  standard 
deviations  are  in  good  agreement  witli  those  for  l62 
(Table  IV)  -  -  an  in-plane  position  error  of  magnitude  40  to 
50  11. m. ,  and  an  out-of-lane  position  error  of  less  than  10  n.m. 
Noting  tliat  a  chord  of  30  n.m.  length  subtends  an  angle  of 
less  than  2**  at  1500  n.m.,  we  see  that  the  orbit  estiiaated 
from  the  initial  conditions  generator  may  be  of  sufficient 
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i-.cc  racy  to  react:uiie  u  satellite  after  one  or  nioro  revolutions 

if  t'le  antenna  is  iiroyramraed  to  remain  fixed  in  inertial  sj-ace 

for  several  minutes.  This  program  might  be  accomplislied  by 

e;:i;ct  co  iputation;  or  a  linear  approxim^ition  may  be  employed. 

La  iaiti\.!.icnt  the  linear  approximation,  a  single  acr,uisition 

Lime  t  ,  .lid  radar  co-ordinates  b  ,  A  ,  E  ,  are  computed, 
a  a'  a  a’ 

llieu  ;;L  time,  t  +  At,  the  radar  is  directed  to  the  co-oidinates 

a 

.3  +  AS,  A  +  AA,  E  +  AE,  .-.here  the  deltas  ui'e  given  by 
a  a  a 

ei.uation  [  3  ]  • 


Itef erence :  Midas  Simulation  and  Evaluation  Study 

Rc^.ort.  iVDL  -  TK1597.  27  ilec.  19ali  Ap]>eniiix  E. 


TABLE  1 

ClLUtiUTEltlSTICs  OF  46l  and  l62  DATA 
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TAiiUi  lA. 
STATION  LOCATIONS 


STATION  NOUTU  LATITUDE  EAST  LONGITUDE 

_ (DEO.  )  _ (DEG  ) 


NEWS 

42.95 

288.37 

VAFB 

54.79 

239.50 

HAWI 

21.57 

201.74 

KODl 

_ iL.60 _ 

207.83 

SKRlAL 

LSTIMArED  FEKIOD 

p  (sec  .  ) 

TABLE  II 

AND  PEIllOD  EHilOHS 

A?  (SEC.) 

FOR  461  RUNS 

P  -  AF  (sec.) 

544 

10519 

+520 

9999 

551 

10002 

+31 

9971 

552 

9862 

-95 

9957 

553 

10054 

+67 

9987 

561 

10112 

+180 

9932 

MEAN 

10109.8 

140.6 

9969.2 

STD. 

DEV.  246.8 

233.7 

26.2 

TABLE  III 

ESTIILVTED  ELhMENTS  FOU  46 1  ItUNS 


Si2ilAL 

A(N.M. ) 

e 

l(U^'.i) 

OiDEG) 

P  (S..C) 

544 

5598.99 

.048837 

95.99089 

265. 612 53 

10519 

551 

541S.93 

.013331 

95. 88709 

265.77391 

10002 

552 

5363.55 

.014089 

95.77398 

265.65816 

9862 

553 

5432.86 

.018370 

95.88710 

265.61156 

10054 

561 

5453. 80 

. 022394 

95-85065 

265-58401 

10112 

ME,\N 

5453.626 

.0-.:34042 

99.877942 

265.647994 

10109.8 

STD.  DEV. 

87.86 

.01468 

.07823 

.07525 

246.8 
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TABLK  IV 


ESTIMATED 

ELEMENTS 

FOR  162 

RONS 

SERIAL 

A  (N.M.) 

a 

1  (DEG. ) 

fi  (deg.  ) 

p  (sec.) 

476 

3736.51 

.01795 

82.40315 

-2.14883 

5734.6 

500 

3793.19 

.02164 

82.42089 

-2.17686 

5865.6 

505 

3702.16 

.02212 

82.28449 

-2.25360 

5655.7 

513 

3740.33 

.01612 

82.24016 

-2.29336 

5743.4 

516 

3768 . 80 

.02265 

82.36774 

-2.22675 

5809.1 

523 

3769.50 

.01613 

82.34294 

-2.28153 

5810.7 

MEAN 

3751.748 

.019435 

82.345228 

-2.230155 

5769.85 

STD.  DEV. 

32.08 

.00305 

_ -0697 

.0577 

_ 21.56 _ 

T/\BLE  V 

FOLK  FA5S  STATISTICb  OF  46l  ELEMENTS 

_ A  (N'  M-) _ e _ I(DEG)  n(DEG.)  ?(SEC)  P-Af(SEC) 

MEAN  5417.285  .017046  99.849705  265.65691  10007.5  9961.75 

STD.  DEV.  38.58  .00420 _ .0733 _ .0838  106.9  23.3 
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